Full Counting Statistics of Charge Transfer in Coulomb Blockade Systems 
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Full counting statistics (FCS) of charge transfer in mesoscopic systems has recently become a 
subject of significant interest, since it proves to reveal an important information about the system 
which can be hardly assessed by other means. While the previous research mostly addressed the FCS 
of non-interacting systems, the present paper deals with the FCS in the limit of strong interaction. 
In this Coulomb blockade limit the electron dynamics is known to be governed by a master equation. 
We develop a general scheme to evaluate the FCS in such case, this being the main result of the 
work presented. We illustrate the scheme, by applying it to concrete systems. For generic case of 
a single resonant level we establish the equivalence of scattering and master equation approach to 
FCS. Further we study a single Coulomb blockade island with two and three leads attached and 
compare the FCS in this case with our recent results concerning an open dot either with two and 
three terminals. We demonstrate that Coulomb interaction suppresses the relative probabilities of 
large current fluctuations. 

PACS numbers; 73.23.-b, 72.70.+m, 05.40.-a, 74.40.+k 



I. INTRODUCTION 

The current fluctuations in the various mesoscopic sys- 
tems have been the subject of both theoretical and exper- 
imental research in the last two decades. Traditionally, 
the attention was focused on the shot noise phenomenon. 
The shot noise is the main fundamental source of current 
noise at low temperatures. In classical systems shot noise 
unambiguously related to the discreteness of the electron 
charge. In quantum system the shot noise can be used as 
unique tool to reveal the information about the electron 
correlations and entanglement of different kind. The in- 
vestigation of the quantum shot noise cross-correlations 
in the multi-terminal mesoscopic devices is the new trend 
in this field which has attracted much attention as well. 
The most achievements in the study of the shot noise 
phenomena have been summarized in the recent review 
articles. 

Alternative way to investigate the current correlations 
in the mesoscopic systems has proposed in the pioneering 
work by Levitov et. ala. This new fascinating theoreti- 
cal approach, known as the full counting statistics, yields 
not only shot noise power but also all possible correla- 
tions and momenta of charge transfer. The essence of 
this method is an evaluation of the probability distri- 
bution function of the numbers of electrons transferred 
to the given terminals during the given period of time. 
The first and the second moments of this distribution 
correspond to the average currents and the shot-noise 
correlations, respectively. The probability distribution 
also contains the fundamental information about large 
current fluctuations in the system. 

Initially, FCS methodi made use of the scattering 
approach to mesoscopic transport. It was assumed 
that the mesoscopic system was completely character- 
ized by its scattering matrix. This method enabled 
the study the statistics of the transport through the 



disordered metallic conductoril and the two-terminal 
chaotic cavityQ. Muzykantskii and Khmelnitskii gener- 
alized the original approach to the case of the normal 
metal/superconducting contacts. The very recent devel- 
opment in this field is the counting statistics of the charge 
pumping in the open quantum dotSErlj'lD. 

The use of multichannel scattering matrix of the sys- 
tem was crucial to obtain the results of the above men- 
tioned works. However, such approach leads to the diffi- 
culties in case of practical layouts, where the scattering 
matrix is random and cumbersome. They become appar- 
ent especially in case of multi-terminal geometry. To cir- 
cumvent these difficulties one evaluates the FCS, with the 
semiclassical Keldysh Green function methodB or with 
its simplification called the circuit theory of mesoscopic 
transportB. The Keldysh method to FCS was first pro- 
posed by one of the authors in order to treat the effects 
of the weak localization corrections onto the FCS in the 
disordered metallic wires. The method proves to be very 
flexible and has been recently aaplied to the FCS in su- 
perconductingJbjeterostructureaij, multi-terminal normal 
metal systemalJ as well as in the three-terminal super- 
conducting beam splitteiO. 

The above research addressed the FCS of non- 
interacting electrons. Since the interaction may bring 
correlations and entanglement of electron states the 
study of FCS of interacting electrons is both challenging 
and interesting. In this paper we present an extensive 
theory of FCS in mesoscopic systems placed in a strong 
Coulomb blockade limit. 

Note, that the shot-noise in the Coulomb blockade cle^ 
vices has attracted the significant attention. KorotkovtJ 
and Hershfield et.aZjlj presented the first theory in the 
framework of "orthodox" approach to single electron 
transport. Later on Korotkov also studied the frequency 
dependence of the shot noise by means of Langevin ap- 
proach both in low (classical) and very high (quantum) 
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frequency limits. The frequency dependence of the shot 
noise in thej-sinffle electron transistor was also investi- 
gated in RefOO The ferromagnetic sip«;le electron tran- 
sistor was considered by Bulka et. al. E3. The-shot noise 
experiments were performed by Birk et. aZ.tj. In this 
work, the nano-particle in between the STM (scanning 
tunneling microscope) tip and the metallic electrode was 
used to form the Coulomb blockade island and the quan- 
titative agreement with the theory of Hershfield et. al. 
was found. 

The electrons dynamics in Coulomb blockade limit is 
fortunately relatively simple. When the cotunncling phe- 
nomena is disregarded, the evolution of the system is gov- 
erned by a master equation. The charge transfer is thus 
a classical stochastic process rather than the quantum 
mechanical one. Nevertheless the FCS is by no means 
trivial and has not been studied yet. In the present pa- 
per we have developed the general approach to FCS in 
the Coulomb blockade regime. This is the central result 
of the paper. Our method turns out to be an elegant 
extension of the usual master equation approach. We 
also predict the FCS in different Coulomb blockade sys- 
tems. The previous results for the zero frequency shot 
noise power can be evaluated in our approach as the sec- 
ond moment of charge transfer probability distribution 
function. 

The paper is organized as follows. In the section II we 
start by presenting the two physical systems to be treated 
within the master equation. Basing on these prototypes 
we formulate the general model. We derive our approach 
to the FCS in the section III. In the section IV we applied 
to the single resonant level and consider its relation to 
the scattering approach. Two- and three-terminal single 
electron transistors are considered in section V. We also 
compare their FCS in the Coulomb blockade limit with 
that one of non-interacting electronsl^. We summarize 
the results in section VI. 



II. SYSTEMS UNDER CONSIDERATION 

The dynamics of many different physical systems can 
be described by master equation. For our purposes it is 
convenient to write it down as follows 

^^\p,t) ^ -L\p{t)) (1) 

where each element of the vector \p{t)) corresponds 
to the probability to find the system in the state n. The 
matrix elements of operator L are given by 

Here r„^,„ stands for the transition rate from the state 
m to the state n. 7„ stands for the total transition rate 
from the state n. Thus defined operator L always has a 
zero eigenvalue, the corresponding eigenvector being the 
stationary solution of the master equation. 
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FIG. 1: The single resonant level system, formed by the two 
tunnel barriers. The resonant level in the quantum well is 
shown by the dashed line. 



A variety of Coulomb blockade mesoscopic systems 
obey Eq. (|^) under appropriate conditions. Here the 
main advantage of the master equation approach is a 
possibility of non-pertubative treatment of the interac- 
tion effects. In what follows, we first remind the master 
equation description of two simple systems: single reso- 
nant level and many-terminal Coulomb blockade island. 
On the basis of these examples we will sketch the master 
equation for the general Coulomb blockade system. This 
will prepare us to the next section where we derive the 
FCS method. 



A. Resonant level model. 

An elaboraljed model of the resonant center was pre- 
sented, in Ref.t3. It was subsequently improved in the 
worked to include the Coulomb interaction. One the 
physical realization is disordM:ed tunneling barrier which 
is placed between two leads. EHl. At sufficiently low tem- 
peratures the main mechanism of transport in this sys- 
tem is the resonant tunneling via localized states formed 
by impurity centers. Another physical realization is the 
Coulomb blockade quantum dot in the low temperature 
regime /cbT ^ AE, where AE is a mean separation be- 
tween the energy levels in the dot.E^ By applying the 
gate voltage, one can tune the given level to be between 
the chemical potentials of the leads. (See Fig.l). We 
consider below two limiting cases where one disregards 
either double occupancy of the level or on-cite Coulomb 
interaction. 

In the strongly interacting case the double occupancy 
of the resonant level is entirely excluded due to the 
Coulomb repulsion U. Then the system can be found 
only in two different microscopic states: one with no elec- 
trons, and another with a single electron. The transport 
through the level can be described by master equation 
approach, provided the applied voltage or the tempera- 
ture are not too low, i.e. max{ey, kgT} ^ ?ir^(/j). Here 
^L{R.) SLie the quantum- mechanical tunneling rates from 
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the left (right) electrode onto the resonant level. We will 
also assume that at the relevant energy scale, given by 
m.eLK{eV, ksT}, the rates ^l{r) are energy independent. 

Under above assumptions the transition rates in Eq. 
(1) are given by 

Ti^o = 2ri/i(e,) + 2rH/7?(e.) (3) 
To^i = ri[l-/i(eO]+rfl[l-/i^(e,)] 

The microscopic states {0} and {1} denote the situation 
with no and one electron, respectively. Fermi function 
fL(R){^) = (1 + exp[(e - ^iL(R))lkT])-^ accounts for the 
filling in the left (right) lead and is the position of the 
resonant level. The factor 2 in the rate Fi^o stems from 
the fact that two quantum states, with spin up and down, 
are available for tunneling. The description in terms of 
rates (^ is correct when the Coulomb repulsion is strong 
enough: U » maxjeV, ksT}. 

The opposite limit is the case of vanishing Coulomb in- 
teraction. In this case the spin up and down channels can 
be treated independently. Each of them can be described 
by the master equation, provided the same condition as 
before is fulfilled: maxjeV, fcsT} ^ KTLiR)- For both 
spin directions the rates are written as 

Ti^o = TLfL{e^)+VRfR{e,) (4) 
Eo^i = ri[l-/i(eO]+rfl[l-/fl(e,)] 

Here the indices {0} and {1} denote the filling factor of 
the level by electon with a chosen spin. 

Disregarding of Coulomb interaction is not adequate 
for a realistic system. However the latter model is worth 
to consider as well. The point is that the statistics of the 
charge transfer in this case can be also evaluated in the 
framework of the non-interacting scattering approachH, 
thus providing the way to establish the consistency of 
two approaches to FCS. 
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FIG. 2: The equivalent circuit of the N-terminal Coulomb 
blockade island. Each junction k is biased by the external 
voltage source Vk. The dot itself is capacitatively coupled 
with the gate voltage Vg. 



B. Many-terminal Coulomb blockade island. 

The electrical circuit incorporating the quantum dot 
with several terminals is shown in Fig. 2. This circuit is 
an extension of the usual sejt-up in case of conventional 
two-terminal quantum dot.E3 At the present stage of nan- 
otechnology the mesoscopic system, associated with this 
circuit, can be realized with the use of two-dimensional 
electron gas (2DEG) in the GaAs/AlGaAs heterostruc- 
tures. 

The essential elements of the circuit shown in Fig. 2 
are the resistances of the contacts, the mutual ca- 
pacitances Cfc between the leads and the island and 
the external dc voltage sources 14. Correspondingly, 
Cg and Vg denote the gate capacitance and gate volt- 
age, which is used to vary the offset change on the is- 
land. We assume that the dot is placed in the Coulomb 
blockade regime, Rk ^ Rq — 27r?i/e^. In order the 
Coulomb blockade effect will be observable the condi- 
tion fcsT ^ Ec — e^/2Cs is also required. Here Ec is a 



charging energy of the island, Ce = Y^^=i^k 
sum capacitance of the system and A'^ > 2 is a number of 
leads attached to the dot. We also assume the tempera- 
ture to be rather high, ksT ^ Ai?, with AE being the 
mean level spacing in the dot, so that the discreteness of 
the energy spectrum in the island is not important. The 
possible effects of co-tunneling will not be discussed in 
the paper. Therefore the characteristic scale of applied 
voltage eV is assumed to be greater than the Coulomb 
blockade threshold, eV > Ec- 

Under the above conditions the multi-terminal dot is 
fairly well described by the "orthodox" Coulomb block- 
ade theory. One can consider the excess number of elec- 
trons on the island (— n) as a good quantum number, 
corresponding to the macroscopic state of the system. 
The tunneling of electrons will occur one by one, increas- 
ing or decreasing the charge Qq — ne on the island by 

(k) 

±e. The corresponding tunneling rate F„_|l]^^„ across the 
junction k is expressed via the electrostatic energy differ- 
ence Ae'^^i^h between the initial (n) and final (n ± 1) 
configurations 



Cg is a 



■ n±l* 



ae. 



(fe) 

n±l<— n 



ik) 

The evaluation of Ai?„j:;^^„ can be dojae along the same 
lines as in the case of two-terminal doto The result reads 



exp[ 



-AE' 



ik) 

n±l< 



-JkBT] 



(5) 



A£;Si^„ = ±e(H-X4)(n) 



2Cj 



(6) 



where Vb(n) is the electrostatic potential on the island. 
It is written as 



, 1 ^ _ 
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Here Ci = Ci + Cg/N and we also assumed that the dot is 
biased such a way, that external voltages are subjected to 
the condition J2iLi — 0- this case the gate voltage 
Vg can be used to influence the offset charge q = CgVg 
on the island in a controlled way. 

Neglecting the quantum correlations between differ- 
ent tunneling processes, we may write down the master 
equation (|l|). It connects the states with different is- 
land charge, the total transition rate from the state n to 
n±l being the sum of tunneling rates over all junctions: 
r _ p(fc) 



C. Universal model. 

We now outline the universal model, which is an ex- 
tension of the preceding two. The physical realization of 
this model can be either an array of Coulomb blockade 
quantum dots or a mesoscopic system with a number of 
resonant levels. We assume that all the relevant condi- 
tions, mentioned previously, are satisfied and therefore 
the description in terms of master equation is valid. It 
is convenient schematically to represent the system as a 
graph (see Fig. 3), so that each node a would correspond 
either to the single dot, the single resonant level or the 
external terminals and the line (a, /?) would be associated 
with the possible transition. Let M be the total number 
of nodes in this graph. We denote nodes by the Greek 
indices and numerate each line by the Latin indices, i.e. 
we will write k = (a,/3), where k = 1 . . .L and L is a 
total number of lines. In case of many-dot system each 
line k unambiguously corresponds to the tunnel junction. 
In case of mesoscopic system with many resonant levels 
one has no a direct one-to-one correspondence between 
the line of the graph and the junction. Still, for the sake 
of convenience, we will sometimes refer below to each 
line of the graph as to the "junction" when it may not 
lead to contradiction. The lines are assumed to be di- 
rected, thus specifying the sign convention for a current 
Ik through the "junction" k. There are N external junc- 



tion k = 1 . . . N, {N < L), connecting the terminals with 
the system. The currents through these junctions are di- 
rectly measurable and hence are of our main interest in 
this paper. 

The macro- or microscopic state of the universal 
model is given by a set of occupation numbers \n) ~ 
\ni, . . . , um)] is equal to any integer for the array of 
quantum dots and refers to the excess charge on the is- 
land a; in case of many resonant levels system Ua = 0, 1, 
since it is assumed that the double occupancy is totally 
suppressed due to the strong Coulomb interaction. 

We now consider the general properties of the master 
equation (|l|) describing the above model. Owing to the 
fact ^„ Lnm — 0, the L operator has the right, |po), and 
the left, (go I, eigenvectors corresponding to zero eigen- 
value 

L|po) = 0, (go|i = (8) 

We assume that they are unique, otherwise it would im- 
ply the situation with two or more non-interacting sub- 
systems. The vector \pq) gives the steady probability dis- 
tribution and (gol = (1, 1: ■ ■ • J !)• Since we are interested, 
in general, only in the permanent, but not the transient 
processes in the system it is naturally to restrict the con- 
sideration only to absorbing states n. Thus we will ex- 
clude all the transient ones n', for which r„„' > but at 
the same time r„/„ — 0. We assume that the L operator, 
bounded to the absorbing states, has a complete set of 
left and write eigenvectors 

L\pk) = bfe)Afe, Xk{qk\ = {qk\L, ^ \pk){qk\ = i 

k 

(9) 

where J is a unitary operator in the absorbing subspace. 
For any physically reasonable L operator Aq = and 
ReAfe > for fc ^ 0. 

For the following it is also useful to represent L oper- 
ator in the form 
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FIG. 3: The graph of universal model (See the main text). 
The terminals are connected with the system via external 
junctions 1,2 and 3. The nodes a, /3 and 7 are either resonant 
levels or dots, linked with each other by internal junctions 
fc's. The arrows denote the conventional direction of a current 
through each junction. 



L 

L = 7-f, f = ^(fW+fi-)) (10) 
fc=i 

where 7 is the diagonal operator in the basis \n) of the 
system configuration and f are associated with the 
tunneling transitions through the "junction" k = {a, 13): 

7 = 5]|n)7(n)(nl, f^^^ =Y.\n')T'^^\n){n\ (11) 
{"} {«} 

The state \n') = |ni, . . . , n^, . . . , n^, . . . , nAf ) results 
from the state \n) by appropriate changing the corre- 
sponding occupation numbers: n'^ = na — <7k, n'^ — 
TT-fj + ak, where ak = ±1 denotes the direction of the 
transition. 
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III. THE FCS IN THE MASTER EQUATION 

In this section we derive the general scheme to evalu- 
ate the FCS in the arbitrary mesoscopic system placed 
in the strong Coulomb blockade regime. In what follows, 
we assume that the mesoscopic system under study can 
be described by the universal model outlined in the pre- 
ceding section. We also assume that the transition rates 
of the corresponding master equation are known. The 
physical origin of these transition rates is a pure quan- 
tum mechanical phenomenon of tunneling. However, the 
reduced description of the system dynamics by means 
of master equation corresponds to the Markov classical 
stochastic process. The mutual quantum correlations be- 
tween subsequent tunneling events are not taking into 
account in such approach. It is a price one has to pay 
to go from the quantum-mechanical to classical descrip- 
tion. On the other hand, it enables us to use classical 
probability methods to derive the FCS. 

In the r&Uowing we will partially use the notation of 
the bookEa. We assume that for each outcome of the 
experiment one can put into correspondence a random 
set of points at the time axis, obeying the condition 



-r/2 > Ti > T2 > • • • > Ts_i > r, > -r/2 



(12) 



Here Ti is the instantaneous moment of z's transition 
in the system, and s = 0,1,... is a non- negative in- 
teger, corresponding to the total number of transitions 
during the time of measurement T. In the end of cal- 
culation this time should be understood as infinitely 
large, T +oo, so that in average s 3> 1. Given this 
set of points we introduce the elementary random event 
= (ti, /ci, (Ti; . . . ; T^, fcs, CTs). It corresponds to the ex- 
perimental outcome, when at time Ti the tunneling hap- 
pens through the junction ki, ai = ±1 being the direction 
of the transition. The events Cs constitute the set ft of 
all possible experimental outcomes. 

At the next step one should define the measure 
(or the probability) d/x(C) at the set ft. For this 
purpose we may very generally introduce the se- 
quence of non-negative probabilities Qs{{Ti,ki,ai}) 
Q(ti, ki, ai; . . . ■,Ts, kg, CTs) > defined at the domain (|lz 
so that 

-t-oo 

(i/^(C) = Qo + ^ E Qs{{n,h,ai\)dTi . . .dTs (13) 

The functions Q are normalized according to the condi- 
tion 



EE /•■■/ Qs{{n.h,<j,})\{dn = i 

(14) 



Each term in Exp. (13) corresponds to the probability of 
an elementary event Cs- 

To complete the preliminaries it is also necessary to de- 
fine the concept of a stochastic process. Mathematically 
speaking, it can be any integrable function A{t) = A{t, C,) 
defined at the space of all experimental outcomes SI and 
parametrically depending on time. It is sometimes con- 
venient to omit the explicit Q dependence. We will use 
a "check" in this case to stress that the quantity in 
question is a random variable. Each stochastic process 
A{t, Q generates the sequence of time dependent func- 
tions ^Ao{t), Ai{t,Ti,ki,ai),..., As{t,{Ti,ki,(Ji})^ Its 
average {A{t))^ over the space Vt is defined as 

{A{t)),,= I A(i,c)dMC) ^^o(OQo + E E 

s=l{ki.a,} 

(15) 

T/2>Ti>--->r,>-T/2 

The analogous prescription should be used, for instance, 
to define the correlations {A{ti)B{t2))Q between any two 
stochastic processes. 

For the subsequent analysis we define the random 
process i'^^\t), corresponding to the classical current 
through the external junction k < N. It is given by 
the sequence of the shot 5-pulses in time 



4=1 

where cr,„ is included to take into account the direction of 
the jump and d{k — ki)= Sk,ki is the Kronecker d symbol. 
Given this defintion at hand, we introduce the generat- 
ing functional S'[{xi(i)}] depending on N counting fields 
Xi(T), each of them associating with a given terminal i: 



exp(-5[{x.(t)}])-(exp{z^ / drx„(T)/(") (r)/e} 



(17) 



The above average is assumed in the sense of defini- 
tions ( p^ ) and (p^. We will refer to S[{xiit)}] as the 
action. Its evaluation is the main goal of this sec- 
tion. The m-order functional derivatives of 5'[{xi(i)}] 
with respect to Xi give the irreducible m-order current 
correlations. First derivatives corresponds to the aver- 
age currents through terminals; the second derivatives 
give the shot noise and noise correlations. In the low- 
frequency limit of current correlations one may use the 
time- independent counting fields Xi- In this case the ac- 
tion 5'[{xi}] allows to express the probability of Ni elec- 
trons to be transferred through the terminal i during the 
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time interval T 



pm}) 



N 

n 



2tt 



-S{{Xi})-^Y.iN,Xi 



(18) 



The above definitions were rather general, but non- 
constructive, since the probabilities Q have not been 
specified so far. One has to relate them with transi- 
tion rates of the master equation. We will show below 
how this can be achieved with use of Markov property 
of the system. Our approach will have much in com- 
mon with the path integral method in statistical physics. 
The advantage of such scheme is that it will readily give 
an explicit algorithm for calculation the generating func- 
tion (0). 

To proceed with the construction of measure ( p^ ) we 
assume that at initial time t = — T/2 the system was in 
the state {n'''^}. Starting from this state it is readily to 
reconstruct the subsequent time evolution of the charge 
configuration {n^^^} *— {n^^^} . . . {n'*"^)} ^ {n^'*)} for 
any given outcome Cs. (See Fig. 3) The choice of (s 
specifies that the transition between neighbouring charge 
states {ri(*~^^} and {n^*)} happens through the junction 
ki — {ai,(3i). Therefore the sequence {n^*)} is given by 
the relation no- — rial ~ ^ki, np, — ny, -\- C/t. , and 

for all 7 7^ Ui and j3i. To write down the 
probability Qs({ri, fci, CTi}) we assume that our event C^s 
constitutes the Markov chain and rely on two facts: (i) 
the conditional probability of the system to survive at 
state n^'-' between the times r^+i and is proportional 
to exp[— 7(n(*))(ri — r^+i)]; (ii) the probability that the 
transition occurs through the junction ki during the time 
interval dri at the moment is given by r[^"'^(n«)dr,. 
These arguments suggest that Q's have the form 

go = Zo"'exp[-7(n(^))T] (19) 
g.({r„fc„a.}) = Zo-iexp[~7(nW)(r/2 - T,)]T[^^'\r^^'^) 
exp[-7(n(i))(Ti - T2)]r(^^)(n(2)) . . . exp[-7(n(^-i)) 
(r,_i - r,)]ri^^)(n(^))exph7(n(^))(r, +r/2)] 

where the constant Zq should be found from the normal- 
ization condition (p^. As we will see below, Zq = 1. In 
Appendix A we show how the usual description of system 
dynamics in terms of master equation can be established 
on the basis of probabilities ([l9|). 

The above correspondence between the random 
Markov chain (s and the probabilities Q's^iol) gives the 
key to evaluate the generating function (|l7). By defini- 
tion (|l^) for any given we have 



measurable current is generated in this case. On averag- 
ing the latter expression over all possible configurations 
with the weight (i/i(C) we see that the generating func- 
tion takes the structure of normalization condition ( |l^ 

Z[{x^{r)}] = cxp(-5[{x.(T)}]) = Qo+ (20) 

S 

Q^{{Ti,ki,ai})Y\_dTi 



E E 



T/2>Ti>--->r^>-T/2 



Here the x-dependent functions Q^{{Ti,ki,ai}) are de- 
fined similar to probabilities (|l^) with the only differ- 
ence that the rates v'j^^n) should be substituted by 

ri^\n) exp{ic7fe XkiTk)} iik<N. 

The expression (|20| ) can be written in the more com- 
pact and elegant way. For that we introduce the x~ 
dependent linear operator defined as 



= 7-rx(^)' 

N 

TxW = E(r 



(21) 



^^+)g«Xfc('r) 



k=l 



in 



/ . k 
k=N+l 



Following the above consideration we have attributed the 
X-dependcnt factor to each operator f (fc = 1 . . . N), 
corresponding to the transition through the external 
junction. The diagonal part and internal transition op- 
erators f with k > N remained unchanged. Then we 
consider the evolution operator U^{ti, associated with 
(^). Since L^{t) is in general time-dependent, U^{ti,t2) 
is given by the time-ordered exponent 

U^{h,t2)^Trexp[- j '(7(t) -f^(r))dr} (22) 

J t2 

The similar construction is widely used in quantum 
statistics. The difference in the present case is that at 
X = the operator U^{ti,t2) gives the evolution of prob- 
ability rather than the amplitude of probability. 

With the use of evolution operator (^) the generating 
function ( pO| ) can be cast into the form 

Z[{Mr)}] = {qo\Ux{T/2, -T/2)\n,) (23) 

To prove it we will argue in the following way. We 
explore the fact, that 7(t) and T{t) commute un- 
der the time-ordered operator in Eq. ( |2^ ) and regard 
T{t) as a perturbation. This gives the matrix element 
{qQ\U^{T/2, -T/2)\ns) in the form of series 

{qo\U^{T/2, ~T/2)\ns) = (go|e"^^|po) + (24) 



expi^^ J dTXn{T)I^''\r,Cs)/e\^l[eM^^^XkM)} E^^ol^- ^^P|- E J - J 



It is assumed here that Xki = if the transition happens 
through internal junction ki > N whereas no physically 



T/2>ri>--->r,.>-T/2 
s 

ri? ^ (ri ) e'^' ^"^^^'^K . .ti^;\Ts) e^^-' ^''^^^'^^\po)l[ dn 
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It follows from the definition ( |l9| ) that each term in 
this series corresponds to the function QJ({Ti, fc^, (7^}), 
namely 



Qa = {qo\e '''^{ris) 



Q: 



K{ri,A;i,cri}) = (go|7Vexp|- /" 7(T)dT| (25) 

J-T/2 ' 



T/2 



Therefore Exp. ( |24| ) and (^3|) are reduced to the pre- 
vious result (2^). This finishes the proof. Note, that 
owing to the property (|^) one has the identity Zq ~ 
(qoI exp(— TI/)|ns) = 1 at x = 0. Therefore the prob- 
abilities ( |l9| ) are correctly normalized. 

The way ( ^3|) the generating function is 
written depends on the initial state \ns) of the sys- 
tem. It looks artificial and we may show that the choice 
of jris) does not affect the final results. It is natu- 
rally to assume that Xfe(t) — > when t — > — oo. i.e. 
physically speaking, the measurement is limited in time. 
To be specific one may suppose that Xfe(t) = when 
-T/2 < t < -T/2 + At. If the time interval At is suffi- 
ciently large as compared with the typical transition time 
then the system will reach the steady state during 
this period of time. It follows directly from the fact that 
exTp{—L At)\ns) — > \po) when At 3> F"^. Thus one can 
substitute \ns) to \pq) in Exp. (p^). Assuming also the 
limit T oo, we arrive to the main result of this section 

+ 00 

exp(-5[{x,(i)}]) = {qo\Trexp[- J L^{T)dT}\Po) 

— OO 

(26) 

We see that the generating function can be written in 
the form of the averaged evolution operator which cor- 
responds to the modified operator L^, depending on the 
counting fields Xi (''")• 

In the rest of the paper we will deal only with the 
low frequency limit of the current correlations, w <C F, 
and concentrate on the particle statistics ([l^). We are 
interested in the probability P{{Ni}) of Ni electrons to 
be transferred through the corresponding terminal during 
a large time interval io S> F"^. In this case one can put 
Xk{t) = Xfe when < t < to and Xk{t) — otherwise. 
The action (Efl) then reduces to the 



S{{x^}) - ioA„,i„({xJ) 



(27) 



where Amindxi}) ^ minimal eigenvalue of the operator 
L^. Thus the problem of statistics in question, provided 
the transition rates in the system are known, is merely a 
problem of the linear algebra. 

In the rest of this section we consider the question of 
the current conservation in the nodes. For this purpose 
we associate the counting fields Xfe with each line k ~ 
{a, /?} of the graph and in the appropriate way modify 
the operator. Then we define the classical current 



operator J^'^' 
the relation 



r{",/3} 



through each line by means of 



= e (Fl+^e''^'= -Dr'e-'^" 



^ dXk 
Its average value can be found via relation 



dS 

to dxk 



4({xa) = f ^ = (9^ |jf b?^) 



(28) 



(29) 



It 



follows from the Eq. 



({xJ) - (gf l^xbD 



(0)\ 



|2 

with {q 



and the fact that 
and being 



the eigenvectors of the operator. 

The average physical currents for k < N are given by 
Ik = Iki{Xi})\^-(f Expanding the vector notation in ( [29| ) 
one gets the usual relation for the current in the master 
equation method. The current is expressed via transition 
rates and the steady probability distribution po{{n}). 

We also introduce the particle number operator h^"^ 
in each node, given by a usual formula 



{»} 

Then after few algebra one see that the relation 



(30) 



(31) 



always holds at any node a. Here the summation is going 
over all nodes /3, connected to a. The choice of the sign 
in front of each term under the sum depends on the situ- 
ation whether the given directed line k = {a, /3} is going 
out or coming into the chosen node a. The Eq. (^ij) gives 
the charge conservation law in the operator language. On 
averaging the latter expression over the steady distribu- 
tion, (gol • ■ • bo), and using (||) and ( p9| ) we arrive at the 
conservation law for the x-dependent currents at each 
node a 



^±/{„,^}({xJ) = 




(32) 



This also ensures the conservation of the physical current 
in the model Ik = 0, where the sum is extended only 
to the external junctions k. It follows from summing up 
the relations (|3^) over all internal nodes a and setting 
X = afterwards. 



IV. RESONANT LEVEL MODEL 

In this section we consider the current statistics of 
the resonant level model. First we focus on the non- 
interacting case. Then we apply the general result of the 
section III to the strongly interacting case and compare 
the statistics in these two different regimes. In the end of 
the section we re-derive the results of the previous works 
concerning the shot noise in these systems. 
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Following the definition (|T|) and the expression for the 
rates (^) , the L^-matrix of the single resonant level model 
in the non-interacting regime reads as 



(33) 



where 



ri^o(x) = TLfLe-''^' + TRfne-'^- (34) 
ro^i(x) = ri(l - /L)e*xi + - /H)e'x^ 

Calculating the minimal eigenvalue of this matrix one 
obtains the current statistics in the form 



S{x)^t„{TL+TR-.JV{x)] 
T^ix) = (rL + rnf + iVLTR X 
[/(_)(e.)(e-^^-l) + /(+)(6.)(e*'^-l)] 



(35) 



Here /(_)(e,) = /L(e,)[l - /i?(e»)], /(+)(e.) = /i^(e^)[l - 
/L(ei)] ^^'^ X — Xi ~ X2- We have also accounted for the 
double occupancy of the level by multiplying the result 
by two. 

Since the Coulomb blockade phenomenon is completely 
disregarded in this model, one might have come to the 
same result in the framework of the pioneering approach 
by Levitov and co-workerso. We will show now that it is 
indeed the case. _ 

Following RcfE the general expression for the current 
statistics through a single contact is written as 



Six) 



-^^/d.ln{l + r„( 



e)x 



(36) 



{fL{e)[l - fR{e)]{e-''^ - 1) + /h(6)[1 - hieW^ - 1)) } 

It is valid for any two-terminal geometry provided the 
region in between two electrodes can be described by the 
one-particle scattering approach and the effects of inter- 
action are of no importance. T'„(e) is a set of transmis- 
sion eigenvalues which are in general energy-dependent. 
Fermi functions include the effects of applied voltage and 
the temperature. In the particular case of a single reso- 
nant level the only resonant transmission eigen-channel 
Trie) plays the essential role, its energy dependence is 
being given by the Breit-Wigner formula 



Trie) = 



(e-Q)2 + (rL + rj^)V4 



(37) 



Here denotes a position of a resonant level. The re- 
sult (36) is more general, than Exp. (^5|). When electrons 



are supposed to be non-interacting the former is valid for 
any temperature. We will show below that one can re- 
produce the statistics ( |35| ) on substituting Trie) into the 
Exp. ( ^6|) and assuming the regime fcsT ^ ?ir. As it was 
discussed previously, this is the condition, when the mas- 
ter equation approach, and hence its consequence (35), 
are valid. 



It is easier to perform the calculation if one first eval- 
uates the x-dependent current /(x) = iie/tQ)dS/dx- It 
reads 



Hx) 



1 



de [/(+)(e)e^>^-/(_)(e)e-*>^] x (38) 



{Tr\e) + [/(_)(e)(e-''^ - 1) + /(+)(e)(e^'^ - 1)] } 

In what follows we assume that the resonant level is 
placed between the chemical potentials /iL{/?,} in the 
leads. Since ksT ^ ^l{r)j the main contribution comes 
from the Lorentz peak and one can put e = in the 
Fermi functions. Therefore we left only with the two 
poles £1(2) — ei ± i\/I?(x)/2 under the integrand (|3^). 
Closing the integration contour in the upper or lower 
half-plane we arrive at 

/(X) = 2eTLrR [/(+)(e,)e^'^ - /(-)(e.)e~'1 / V^) 



On integrating it over x one finds for the Six) — 
ito/ie) Jq Iix')dx' the result (35) obtained by means of 
master equation. Thus, we have verified the correspon- 
dence between two approaches to statistics in the non- 
interacting regime. 

To proceed we address the strongly interacting regime. 
In this case the (2x2) L^-matrix is formed with the use of 
rates (^. It has a structure similar to Eq. (^4[), provided 
the x-dependent rates are written as 



ri^o(x)-2rL/Le-^>^i 
ro^i(x) = ri(i-/i)e*xi 



2rRfRe-'^- (39) 



On evaluating the corresponding eigenvalue one can 
write down the expression for statistics in the strongly 
interacting limit 

Six) = (io/2){ri[l + /i(6.)] + Tr[1 + fRie.O] - 

Vix) = {Tl[1 + fLie^)] + Tr[1 + /«(£,)]}' + 8TlTr x 
[/(_)(eO(e^'^-l) + /(+)(60(e-^>^-l)] (40) 

To proceed we consider the shot noise regime eV ^ 
ksT and assume that the voltage is applied to the right 
electrode as shown in Fig. 2. Then the temperature fluc- 
tuations become non-essential and both statistics ( |35| ) 
and (^) take a rather simple form 

Six)\u=a = ta^TL + Tr- 

V(rfl-rL)2 + 4rirfle-'x| (41) 
ito{2rL + rfl- 

^i2TL-TR)^ + 8TLTRe-^x^ (42) 



^(x)l 



Given the latter expressions at hand one can eas- 
ily re-derive the known result for the average current 
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FIG. 4: The current statistics through the single resonant 
level. 1 and 4 - non-interacting model; 2 - interacting model; 
4 - Poisson statistics 



V. THE FCS IN THE COULOMB BLOCKADE 
QUANTUM DOTS. 

In this section we discuss the application of the method 
to the many-terminal Coulomb blockade quantum dot. 
The consideration will be limited to the two- and three- 
terminal layouts. Our treatment will be mainly numer- 
ical, though some analytical results in the two-terminal 
setup are also plausible. In the beginning of the sec- 
tion a few technical details, which are common for both 
cases, are given. In particular, we establish the relation of 
the FCS approach with the preceding papers, concerning 
the shot noise in the conventional two-terminal Coulomb 
blockade dots. In the following we consider the FCS, 
first for two-terminal, and then for three-terminal config- 
urations. We will also compare the FCS in the strongly 
Coulomb blockade limit with our recent results, concern- 
ing the FCS in many-terminal chaotic quantum dot with 
contacts being the tunnel junctions. 

A. General remarks 



/ — {ie/to)dS/dx\x=o and the shot noise power S'shot = 
{e'^ /tQ)d'^ S / dx'^\x=o in these models. It is conventional 
to represent S'shot in the form S'shot = elF, where F 
is the so-called Fano factor. Then one obtains F = 
(r| -I- r|j.)/(rL -l- Fa)^ in the non-interacting regime and 



F = 
limi 



4r|i,)/(ri + 2rfl)^ in the Coulomb blockade 



As one can see from the Eq. (|4l|), at low temperatures, 
the difference of statistics in the large U limit from the 
one in the non-interacting case is an effective suppression 
of Tji rate by a factor of two. To find the probability 
distribution P{N,to) one can estimate the integral ( p^ ) 
by means of steepest descent method. It is applicable 
in the given case of low frequency regime uj T, which 
we consider, since both the action S{x) 3> 1 and the 
average number of transmitted electrons N — Ito/e ^ 1. 
Then one has to find the saddle point x* of the function 
f2(x) — S(x) + ixIto/e, which is defined by the equation 
/ — {ie/to)[dS/dx]- It is turned out that x* always lies 
on the imaginary axis. This equation can be regarded 
as a parametric relation between / and x* j and with the 
exponential accuracy we obtain the estimation for the 
probability P{T) ^ exp[— Jl(x*)]. 

The results for the statistics (^l]) are shown in Fig. 4. 
The statistics are compared with the Poisson type statis- 
tics S{x) — 2rio[exp(ix) — 1] with the effective rate given 
by F"! = F^^ + Fj^\ Both statistics (|l]) approach the 
Poisson one, provided the system is strongly asymmetric. 



In case of Coulomb blockade island the macroscopic 
state of the system is characterized by the excess charge 
Q = ne, which is quantized in terms of electron charge 
(— e). As we have pointed out in the section II, within 
the " orthodox theory" , the charge Q can be changed only 
by ±e in course of one tunneling event. Therefore the 
master equation connects the given macroscopic state n 
only with the neighboring states nil. The corresponding 
rates r„±i^„ of these transitions are equal to the sum 

(k) 

of iV = 2 or 3 independent probabilities ^n±i^n through 
the different junctions, those are given by (o) and m). 
Along with the lines of section III, in order to find the 
FCS of the charge transfer through the island, we have 
to modify the rates F„-|-i^„ into x-dependent quantities 
^n±i*-n in accordance with the rule ( pT| ) and to evaluate 
the minimal eigenvalue Amin of the three-diagonal matrix 
afterwards. In the given case it is convenient to write 
down the latter problem as the eigenvalue problem for 
the following set of linear equations 



(A - 7„)p„ + rj^^„+ip„+i + F^^„_ip„_i = (43) 
where 7„ = r„^„_i + r„^„+i, and r^^„_,_;^ = 

J2k=i riii^n^^*'^'' ■ Sign index ± denotes the outcome 
(income) of an electron from (to) the island. 



In general, we have treated the problem (43) numer- 



ically. At sufficiently low temperatures fcsT <C i?c, 
which is mainly the case of the following discussion, 
the temperature dependence in rates (|^) is non-essential. 

Then one can set F^^^) = AEl%^^/{e^Rk) when 



AEl%^^/{e'^Rk) > and T[^\n) = otherwise. Thus 
defined rates are linear functions in n. The possible set 
of {n}, corresponding to nonvanishing rates, is limited 
to some interval rimin < n < Umax- Hence Eq. ( ^3| ) be- 
comes a finite linear problem. At higher temperatures 
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ksT < Ec we have found that the increase (decrease) of 
both riinin and n^ax to extra 7-^8 states gives the results 
up to 10~^^ degree of accuracy in course of the numerical 
procedure. 

The matrix of Eq. ( p3| ) is non-Hermitian. This fact 
may cause an instability in the numerical algorithm when 
the range [rimin , ?^max] is large. However, in most practi- 
cal cases this problem can be circumvented by transform- 
ing to Hermitian form. First we note, that one only 
needs to work with pure imaginary counting fields Xki 
as long as the probability P{{Ni} ,1^) is estimated in the 
saddle point approximation. (See the discussion in the 
end of the section III.) Hence the rates r^^„_|_j^ become 
the positive real numbers. Then we can apply the linear 
transformation p'^ — AnPn- It leads to the rates in the 



-V/2 




FIG. 5: The equivalent circuit of the two-terminal Coulomb 
blockade island. 



new gauge F'^^„±i = A„+iF^^„±i^„^ The unknown ^bT <C eVCs the x-dependent rates F^^„±i reads as 



AnS may be chosen in a way that the symmetry relation 
T'^^j^^i = would hold. This gives the recur- 

rent relation An+i/An = ir^<_„+i/F^+i^„)i/2. with 
the use of latter the Eq. (|4^) takes the Hermitian form 
when it is written in terms of p'^ and transformed rates 



-nil 



y^"^. The diagonal term 7„ 



is not affected under this transformation. 

Let us also discuss a useful relation for the shot noise 



correlations Skm = {e'^/to)d^S/dxkdX' 



x=o- 

use the identity Ainin({Xi}) = ilx' \^x\Px'') ™ order to 
express them in terms of eigenvectors , | and eigen- 
values A„ of the matrix L. With the use of standard al- 
gebra and assuming the normalization {qn\Pn') = Sn,n' 
we can cast Skm in the form 

Skm = e^{qo\S^'''"'^\po)+ (44) 
e^X! T" {(lo\Jk\pi){qi\Ji\po) + {qo\Jm\pi){qi\Jk\po)^ 



One can 



l>0 



where the Jk operator was defined by 
d^LJdxkdx 



and S'C^^") = 
holds in 



^^Q. Note, that the relation (|44, 

any basis. One may, in particular, use it for the basis p'„ 
discussed above. In this case one must define the matrix 
elements of Jk as [Jfc]„_^^ „ = {id / dxk)^'n^n±i and to 
evaluate S''-'^'™-' in the same manner. 

The relation (^) represents the fact, that the shot 
noise correlations are defined by the whole spectrum of 
the relaxation times t^^ = Xk in the system. In the 
case of two-teijtapal geometry it coincides with preceding 
results of Ref.BEl. 



B. Two-terminal Coulomb blockade island. 

The electrical circuit with a two-terminal Coulomb 
blockade island is shown in Fig. 5. The dot is biased 
such a way, that V2 = —Vi = V/2. At low temperatures 



pX 
n+li-n 



pX 

n — l<—n 



CiV 

e 

C2V 



C9V9 



R2 Cs 



«Xl 



■(45) 



where Ci(2) = ^i{2) + Cg/2 are effective capacitances. 
The gate voltage Vg can be used to control the offset 
charge go = CgVg on the dot. It can be varied con- 
tinuously in the range — e/2 < qa < e/2. The result- 
ing dimension of the matrix is given by the number 
of absorbing states nmax — ?^min, where nmin (jT-max) is 
the maximal (minimal) integer closest to the points rti^2 
where the rates ^^^i^n vanish. 

First we briefly consider ±he voltage dependence of the 
shot noise in the systenJiJ'El. It was calculated with the 



o 
c 

to 




-1.0 



-0.5 0.0 0.5 

Offset charge (q ^e) 



FIG. 6: The Fano factor versus the offset charge in the 2 - 
terminal Coulomb blockade dot. The parameters are: Ci = 
C2, Ri = IQR2, T = O.OleVsCE; (1) - F = l.Se/Cs, (2) - 
F = 2.0e/CE, (3) - F = 4.0e/CE, (4) - F = 6.0e/CE. 
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use of Exp. (44). The results for the noise-to-current 
ratio (Fano factor) are presented in Figs. 6-8. The 
Coulomb blockade features are strongly pronounced in 
case of asymmetric junctions only. In the experimen- 
tal situation, when the dot is made up with the use of 
2D-elcctron gas in the semiconducting heterostructure, 
the resistances of the contacts are much easier to vary 
than the mutual capacitances. Therefore we have cho- 
sen Ci = C2 and plotted the Fano factor for different 
values of ratio R2/R1 and offset charge go- The curves 
are truncated below the Coulomb blockade threshold, 
where the considered "orthodox" theory is not applica- 
ble. At high values of the ratio R2/R1 they exibit the 
strong characteristic Coulomb blockade oscillations. The 
special points at the voltage dependences occur when 
either nmin or rimax are changed by 1. At high bias 
voltages the noise-to-current ratio saturates to the value 
F = {R\ + Rl)/{Ri + R2Y independently of the capac- 
itances Cfc and the offset charge go- (See the discussion 
below as well.) An increase of a temperature leads to the 
smearing of oscillations due to the additional thermal 
noise. The above results coincide with those, obtained 
previously by Hershfiled et. a/.Ej 

To proceed we turn to the question of the FCS. For 
the sake of clarity we first present our recent analytical 
results for the FCS of the chaotic quantum dot with two 
tunnel junctions, when their resistances Rk <C irfi/e^. In 
this case the effects of Coulomb interaction are negligible. 
Then we trace the differences in the FCS, when the dot 
is placed in the strongly interacting regime Rk 3> TTh/e^ . 

In the non-interacting limit the action S{x) is pt,- 
pressed via the voltage V and the resistances Rk onljO. 
At low temperatures ksT <C eV it has a form similar 
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V/fe/C^] 



FIG. 7: The Fano factor versus the applied voltage in the 
2-terminal Coulomb blockade dot. The parameters are: Ci = 
C2, ksT < eVsCs, go = 0; (1) - R1/R2 = 1, (2) - i?i/i?2 = 
2, (3) - R1/R2 = 4, (4) - R1/R2 = 10 
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FIG. 8: The Fano factor versus the applied voltage at differ- 
ent temperatures. Parameters are shown on the plot. Tem- 
perature is given in terms of charging energy. 



to(gl|) 



2e 



Rn 



(46) 



- i?^^)2 + 4(i?ii?2)-ie»te-xi) I 



It would be completely equivalent to the statistics of the 
charge transfer by non-interacting particles through the 
resonant level if one regards the ratios F^.^ — V/{eRi^2) 
as the effective tunneling rates. The generating func- 
tion (46) gives the above mentioned value F = {Rf + 
Rl)/{Ri + R2f for the Fano factor. 

In the strongly Coulomb blockade limit the action 
5(x), in general, remarkably deviates from Exp. (|4^). 
Still, there are two exceptions, when S{x) resembles the 
statistics (|l]) and (||). 

The first case occurs at low voltages, slightly above the 
Coulomb blockade threshold value, when only one charg- 
ing state is available for tunneling. This situation can 
be easily realized in the asymmetric dot with R2 ^ Ri- 
Then mere two states with n = and n = 1 are involved 
and is reduced to the 2x2 matrix (33). The only 
difference is that the rates F^ contain the voltage de- 
pendence as given by Exp. (ll). Thus the action S{x) 
reproduces the result (pl|), where the rates Tj^ j^ are as- 
sumed to be voltage dependent. 

To proceed we describe the second exceptional situa- 
tion when the action S(y) can be found analytically. Let 
rii_2 be zeros of rates (p5|), i.e. F^_j_;^^„ (711^2) = 0. We 
now interested in the situation when both zeroes ni_2 si- 
multaneously become integers. This situation may occur 
at the limited number of special points Vk at the Coulomb 
blockade staircase when the ratio C1/C2 is close to a ra- 
tional along with the special choice of the offset charge qq. 
(E.g., for the configuration Ci = C2, shown in Figs. 6, 8, 
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this is the case when (i) qq = 0, Vk = {2k + l)e/Cs and 
(ii) Qo = ±e/2, Vfc = 2k{e/C^) with k being integer.) In 
this situation we may show analytically (see Appendix B 
for the proof), that the action S{x) at points Vk takes 
the form similar to the statistics in the non-interacting 
regime (^). The only difference that the voltage V has 
to be substituted to {\Vk \ — e/Cs). The reduction of V 
by the amount of the threshold voltage value e/Cs is 
thus the manifestation of the Coulomb interaction. The 
given statistics is also valid as a limit at high voltages 
V ^ e/Cs. One may conclude it from the physically 
reasonable arguments that the result for the action in 
this limit should be linear function in voltage and must 
not depend on the capacitance ratio C\/C2- Hence, the 
statistics is insensitive to the fact, whether ni.2 are inte- 
gers or not. This also explains the saturation of the Fano 
factor at Figs. 6-8 to the non-interacting current-to-noise 
ratio. 

To access the general situation, the whole problem has 
been treated numerically. We have evaluated the proba- 
bility ( |l8| ) in the saddle point approximation along with 
the same lines as it was done for the resonant level model. 
The Eq. (|9|) was used as the parametric relation between 
the current / and the counting field X = X2 ~ Xi ™- the 
saddle point x*- In Fig- 9 we give the example of the 
logarithm of probability distribution P{I) for the num- 
ber of different voltages V and offset charge go = 0. All 
curves are normalized by the reduced value of the voltage 
{V — e/Cs). For the value V = Se/Cs the renormalized 
logarithm of probability coincides with the one, obtained 
from the non- interacting limit (p6|). We have also plot- 
ted the same statistics in Fig. 4 in the resonant level 
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FIG. 9: The statistics of current in the two-terminal quan- 
tum dot. Parameters are shown on the plot. Curves 3 and 5 
coincide but correspond to different axes. 1) - = 1.5e/CE, 
(2) - V = 2.0e/CE, (3) - V = 3.0e/CE, (4) - = 4.0e/CE, 
(5) - non-interacting regime. 




FIG. 10: The equivalent circuit of the three-terminal Coulomb 
blockade island. 



model for the ratio of rates F^/F/j = 10, when the in- 
teraction effects are disregarded. We see, that in general 
the probability distribution is strongly affected by the 
Coulomb blockade phenomenon, as compared to the non- 
interacting regime. It approaches to the non-interacting 
limit only at rather high voltages V ^ e/Cs. 



C. Three-terminal Coulomb blockade island 

The electrical circuit with a three-terminal Coulomb 
blockade island is presented in Fig. 10. It is biased by 
three external voltage sources so that the current, flowing 
through the third terminal, would split into the first and 
the second ones. The voltages J7i(2) are used to control 
the bias between the 3d and the 1st (the 2nd) terminals: 
Ui{2) = — Vi(2) . The third terminal is biased at voltage 
V3 = (Ui + C/2)/3 with respect to the ground. Such type 
of setup assures the condition Vi + V2 + V3 ~ 0. Hence 
the gate voltage, Vg, can be used as before to control 
the offset charge go — CgVg on the island. As in the 
previous subsection we discuss the only low temperature 
regime ksT <^ e^/Cs. In what follows it is assumed that 
C/2 > C/i. Then, according to Exp. (|-0) and (|2l|) the X" 
dependent rates of the system are written as follows 



F^ 

n — n 



F^+^ (n)e^» + F^+^ (n) 6l(m - 1/2 - n)e'^^ 



in)e- 



1X2 



(n) e{n - 



where 



Tf\n) 



CgVg 



± 



Rk Ce 



1/2)6-'^' 
(47) 



(48) 



and 



,(+) - 



CiUi + C2U2 



(C'i + C3)C/2-CiC/i 



± 



C2U2 - {C3 + C2)Ul 
ei?iCE 
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The point 
1/2) 



is determined by the relation t[ \m 



1/2) — (m is non-integer in general). 
The dimension of the i^-matrix is equal to rimax — 'T-min, 
where Hmax('T-min) Can be found from the conditions 



.(-) 



in) 



> 



{T[+\n) 



^ 0)- 

Ck are defined as Cfe = Cfc -f 
We can see from the Exp. 



The effective capacitances 
C,/3 

(p7|) that there are four ele- 



mentary processes of charge transfer in the system at low 
temperatures, each of them is being associated with the 
pre-factor e^*'^'' . The presence of the exponents e'^^ and 
e~^^^ corresponds to the charge transfer from the third 
terminal into the island and from the island into the sec- 
ond terminal, respectively. Hence, the random current 
through the 3d (2nd) junctions always have the positive 
(negative) sign. Two factors e^^'^ stems from the charge 
transfer through the first junction in the direction either 
from the island into the first contact or vice versa. There- 
fore the random current Ii is able to fluctuate in both 
directions. 

Let us first consider the shot noise correlations in the 
system. For that it is useful to introduce the (3 x 3) ma- 
trix F with elements Fkm — Skm/sI'E, where the current 
correlations Skm are given by (Q) and /e = 
The matrix f is a generalization of the Fano factor for 
the multiterminal system. It is symmetric and obeys 
the relation X]i=i = 0. It follows from the general 
law of the current conservation in the system. For the 
considered 3-terminal dot we have found that the cross- 
correlations Fkm {k ^ m) are always negative for any set 
of parameters. 

In Fig. 11 we give the illustrative example of the volt- 
age dependence of the shot noise correlations Fkm for the 
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FIG. 11: The matrix F of auto- and cross- shot noise correla- 
tion versus voltage Ui for the 3-terminal quantum dot setup. 
Parameters are shown on the plot. (1) - (2) - I-F12I, (3) - 
|fi3|, (4) - F22, (5) - IF23I, (6) - Fs3. 
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FIG. 12: The matrix F of auto- and cross- shot noise cor- 
relation versus the offset charge for the 3-terminal quantum 
dot setup. Parameters are the same as on the Fig. 11. The 
voltage [/i = C/2/4 = 1.25e/CE. (1) - i^ii, (2) - \F^-,\, (3) - 
iFial, (4) - F22, (5) - IF23I, (6) - 



certain choice of parameters. (For the cross-correlations 
the modulus |i^fcm| are given.) As in the 2-terminal case 
the Coulomb blockade features are strongly pronounced 
only for the asymmetric setup. The results in Fig. 10 
corresponds to Ri = = R2/W, C\ = C2 = C'a and 
U2/U1 — 4. The latter ratio of voltages has been cho- 
sen on the basis of arguments that for a given value of 
resistances Rk it would split the average current 1^ into 
two equal currents = I/2I = h/2 provided one could 
apply the usual linear Kirchgoff rules to this circuit. In 
Fig. 12 the dependence of the shot noise correlations on 
the offset charge is shown for the same set of parameters 
and the value of Ui = 1.25 e/Cs- The special points of 
both these dependences occur when either n„iin) '^max or 
the integer part of m are changed by ±1. As the result 
we observe multi-periodic Coulomb blockade oscillations 
in the offset charge dependences in contrast to the single 
periodic oscillations in the two-terminal case. 

We now proceed with the consideration of the FCS. As 
before, the action S{{xi}) has been calculated with the 
use of (27) and afterwards the probability ( |l8| ) has been 
estimated by means of the steepest descent method. The 
difference with the two-terminal geometry is that three 
currents Ik and three counting fields Xk are now involved. 
Due to the current conservation Ik — Q for any 
plausible fluctuation, only two currents are independent 
and thus the action S{{xi]) depends on the differences 
Xi] — Xi~ Xj only. In what follows we have chosen Ii and 
I2 as the independent variables to plot the logarithm of 
probability lnP(Ji,/2). With the exponential accuracy 
it is given by lnP(/i,/2) ^ e^^''^ \ where x* is a saddle 
point of the function fl{x) = S{x)+iXihtQ/e+iX2l2to/G- 
The results for lnP(/i,/2) are shown in Figs. 13, 14, and 
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FIG. 13: The logarithm of current distribution In P{Ii, I2) 
in the 3-terminal quantum dot as a function of current Ji, 
under condition I2 = {^2)- Parameters are shown on the 
plot. (1) - f/i = 1.25 e/Cs, (2) - Ui = 2.0e/Cs, (3) - Ui = 
4.0e/CE, (4) - Ui — 10.0 c/Ce; curve (5) corresponds to the 
non-interacting regime. 



15. From the contour map on Fig. 15 we see that P(/i, 12) 
is non-zero in the quarter Ii < 0, I2 < of a current 
plain (/i, I2) and in the region Ii < I/2I belonging to the 
quarter Ii > 0, I2 < 0. This range of plausible current 
fluctuations results from the x-dependence of rates ( ^ ) 
and the restriction Ik = 0. As it was discussed above, 
any current fluctuation obeys the condition I2 < and 
I3 > 0. There is no additional restriction to the latter 
when Ii < 0. On the other hand, it follows from the cur- 
rent conservation, that Ii = I/2I — ^3 < 1^2 1 in the region 
h > 0. 

It is also worth to compare the above results with the 
FCS in the 3-terminal chaotic quantum dot when its con- 
tact are tunnel junctions with resistances R'j^^ ^ e^/irh. 
In this limit the effects of interaction are negligible and 
electrons are scattered independently at different ener- 
gies. Provided U2 > Ui the generating function S{{xi}) 
in the given case is a sum of the two independent pro- 
cesses 

'S'(Xl,X2,X3) = <S'l(Xl,X2,X3) + S'2(Xl,X2,X3) (49) 

Here 

'5'i(xi,X2,X3) = ~^{^^ +G2 + G3 - 
^/{Gl + G2 - G3)2 + 4G3e'X3(Gie-«xi + G2e-«X2)| 

52(xi,X2,X3) = ^£i:^{g, +G2 + G,- 
\/(Gi -FGa -G2)2 -h4G2e-^X2(Gie^xi +036^x3)^ 



and Gk = Rf. are the conductances of the junctions. 

The logarithm of probability lnPo(-^i5 -^2), evaluated 
with the use of statistics (^9|), is shown by the dashed 
line in addition to the previous curves in Figs. 13, 14. 
Its contour map for the same values of parameters is 
also separately presented in Fig. 16. The maximum of 
lnPo(-^i> -^2), as expected, occur at Ii = I2 = Ui/3Ri. 

The two main conclusion which one can derive on com- 
paring the statistics in the two limiting regimes (strongly 
interacting in case of the Coulomb blockade quantum dot 
and the non-interacting one in the open chaotic quantum 
dot) are as follows. First, we see, that in spite of the 
different regimes, the qualitative dependence of proba- 
bilities versus the currents is similar for both statistics. 
Second, we may conclude, that the relative probability 
of the big current fluctuations in the Coulomb blockade 
limit is suppressed with respect to the situation in the 
non-interacting regime. This is a characteristic feature 
of both two- and three-terminal quantum dots. This be- 
havior is rather natural. It stems from the fact, that any 
big current fluctuation in Coulomb blockade dot is re- 
lated with the large accumulation (or depletion) of the 
charge on the island. The latter costs the extra elec- 
trostatic energy and hence the probability of such fluc- 
tuation is drastically decreased. The suppression of the 
average current due to the presence of the Coulomb gap 
in the interacting system is the particular example of this 
behavior. 




FIG. 14: The logarithm of current distribution \n P{Ii, I2) 
in the 3-terminal quantum dot as a function of current I2, 
under condition Ji = (Ii). Parameters are the same as in 
Fig. 13. (1) - Ui = 1.25 e/CE, (2) - f/i = 2.0e/CE, (3) - 
Ui = 4.0e/CE, curve (4) corresponds to the non-interacting 
regime. 
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FIG. 15: The contour maps of the current distribution 
log[P(/i, /2)] in the S-terminal Coulomb blockade dot. Pa- 
rameters are the same as in Fig. 13. Ui — f/2/4.0 = 1.25 e/Cs- 
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FIG. 16: The contour maps of the current distribution 
log[P(Ji, 12)] in the 3-terminal chaotic quantum dot with tun- 
nel contacts. Ui = (72/4.0 = 1.25 e/Cs- 



VI. CONCLUSIONS 

To conclude, in the present paper we have developed 
the constructive scheme to evaluate the FCS of charge 
transfer in the Coulomb blockade systems. This scheme 
is rather general and universal and is applicable to any 
strongly interacting system, provided the latter can be 
described classically in the framework of the master equa- 
tion approach. The method proposed consists in the 
transformation of the initial linear operator L of the mas- 
ter equation into the auxiliary ^-dependent linear oper- 
ator L^. Each non-diagonal term of this new operator, 
associated with the particular transition in the system, is 
modified by the exponential prefactor e^*-^'' in order to 
take into account the electron jump through the junction 
k during the tunneling event. The generating function of 
the charge transfer through the whole system is then pro- 
portional to the minimal eigenvalue of the operator L^. 

We have applied this scheme to study the FCS in two 
different systems. For a generic case of a single resonant 
level model we have established the equivalence of the 



new method with the scattering approach to the FCS, 
when the particles in the system are non-interacting. 

Afterwards we have considered the FCS and the 
shot noise correlations in the two- and three-terminal 
Coulomb blockade quantum dots. The consideration was 
limited to the temperature regime when the orthodox 
theory of Coulomb blockade phenomenon is applicable. 
For the case of two-terminal dot we have re-established 
all the known results for the shot noise in this system. In 
the three-terminal case we have shown that the auto- and 
cross- shot noise correlations exhibit the characteristic 
Coulomb blockade oscillations as the functions of the ap- 
plied voltages and the offset charge. We have considered 
the question of the FCS as well. In general situation we 
evaluated the probability distribution numerically. How- 
ever, at some special values of parameters in two-terminal 
dot we have managed to find FCS analytically. In these 
exceptional cases the FCS resembles the statistics of the 
charge transfer through the single resonant level. Then 
we compared the statistics in the Coulomb blockade dots 
with our previous results concerning an open dot with 
two and three terminals. We found that the Coulomb 
interaction suppresses the relative probability of the big 
current fluctuations. 
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APPENDIX A 

In this Appendix we show that the construction of 
the probability measure on the basis of Markov chains 
Cs, which was used to derive the main result of section 
III, leads to the usual description of the system dynam- 
ics in terms of master equation. This correspondence is 
achieved in the standard way of probability theory by 
introducing the stochastic process h(t) corresponding to 
the island charge at a given time t 

s 

n(t,C.) =n,, +^CT,0(<-T,) (Al) 

i=l 

Similarly one can consider the random number of elec- 
trons n^'''>{t) transferred through the junction k after 

t > -T/2 

S 

i=l 



16 



The random variables h{t) and fi'^^\t) are subjected to 
the relations 



N 

k=l 



(A3) 



After that we can introduce the probability distri- 
bution P{n, t) and the joint probability distribution 
P(ni, ii; 7^2, ^2) (^1 > h) of the process n{t) 



P{n,t) 



(A4) 



P{ni,ti;n2,t2) ^ / S{ni - n{ti, C))S{n2 ~ n{t2, C))dKC) 
Jn 

Their ratio P(rii, ti| ^2, ^2) = P{ni,ti;n2,t2) / Pin2,t2) 
gives the conditional probability to find the system at 
state ni at time ti , gi ven that at time t2 it was at state 
712. The integrals (A4) can be efficiently evaluated along 
with the same reasoning as we have used to prove the 
normalization condition. As the result one ends up with 



P(ni,ti| ^2,^2) = {ni\U{ti,t2)\n2) 



(A5) 



The latter expression is the usual way to describe the 
system in terms of master equation. The conditional 
probability P(ni, | n2, ^2) regarded as a function of ni 
and ti obeys this equation with the initial condition 

P{ni,ti) = 5„i,n2 at ti = t2. 



APPENDIX B 

This appendix contains the derivation of the action 
S{x) at low temperatures at some special points 14 in the 
Coulomb blockade staircase in the two-terminal dot. We 
introduce the notation T{^\n) — rj^_|_^„, that enables 
to write down Eq. ( ^ ) in the form 

(A - 7„K + r(-)(n + iK+i + r[+\n - lK_i - 



(Bl) 



If all Xk = then the stationary solution of this equation, 
corresponding to A = 0, satisfies the detailed balance 
condition pn+iTQ ■'(n + 1) = pnr''^\n). In gen eral situ- 



ation, when Xk 7^ 0, one may try to resolve (Bl) making 
use of the substitution 



Pn+i ^ r^'^V) Pn^i ^ y'^i \n) 

Pn yr'^\n + iy Pn r^+^(n-l) 



(B2) 



with unknown constant y to be found. This reduces the 
difference equation (Bl) to the relation 



(A - 7„) + {n)y + r^+) (n)y-i = 0. (B3) 

Here 7„ — Fg (n) -t-FQ^' (n), and F^^' (n) are linear func- 
tions in n, given by Exp. (45). Then one might find the 
two unknown y and A^ on comparing the constant and 



linear in n terms in relation (B3). It yields 



Mx) = Yjy - e/cs) (i?r' + - V^) (B4) 



It looks like we have found in such a way the required 
solution. However, it is not valid at all possible val- 
ues of parameters. The matter is that in the most 
general situation the expressions (^) are not correct 
at points n = rimin and n = nmax. One have to set 



.(-) 



(^max) 



: nmin and n 

Tx (^min) — 



by hand that breaks the 



analytical n - dependence of equations (Bl) and (B3) at 
the boundaries. The only exceptional situation, when 
the above way to solve the Eq. (Bl) is indeed true, corre- 
sponds to the case ni = n^i^, rt2 — ^imax, with 711(2) be- 
ing the zeros of functions F^*' (n). In this case the substi- 
tution (B2) gives p(ri„iax + l) = p(7imin — 1) — and hence 



the actual values of F^ '(^imax + 1) and T)^' {riaiin — 1) 
in Eq. (Bl) play no role. Then we arrive to the action 
S{x) — ioA(x) in the form which was claimed in section 
V(a). 
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